Effect of band structure on field emission of crystalline graphite 
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Abstract 

The field emission of crystalline AAA graphite is studied within a simple analytical approach with 
account of the exact dispersion relation near the Fermi level. The emission current is calculated 
for two crystal orientations with respect to the applied electric field. It is found that the exponent 
of the Fowler-Nordheim equation remains the same while the preexponential factor is markedly 
modified. For both field directions, the linear field dependence is found in weak fields and the 
standard quadratic Fowler-Nordheim behavior takes place in strong fields. A strong dependence 
of the emission current from the interlayer distance is observed. As an illustration of the method 
the known case of a single-walled carbon nanotube is considered. 
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I. INTRODUCTION 



Different carbon-based structures are considered as promising electrode material for field 
emission (FE) cathodes. In particular, the field emission properties of single-walled (SWNT) 
and multi-walled (MWNT) carbon nanotubes [l| as well as graphite films {2} are presently 
under intensive experimental and theoretical investigations. In experiment, many factors 
such as inhomogeneities at the cathode surface, surface contamination (surface adsorbates 
and oxides), local electric fields and barriers, electronic structure of cathode, etc. can dras- 

n 

tically change FE results |3J . In addition, these factors vary from one experiment to another 
thus markedly complicating the theoretical description. Nevertheless, the electronic charac- 
teristics of cathodes should be equally manifested in different experiments. For this reason, 
the effect of electronic structure on the emission features of cathodes is of definite interest. 
For SWNTs this problem was studied numerically in 

BAB by using an approach which 

can be called SIS db method of independent channels. 

In this paper, we present a rather simple modification of this method to study analytically 
the influence of the 3D band structure on the field emission current (FEC) of crystalline 
graphite. As an illustration, we consider the case of SWNT where our approach allows us to 
reproduce the FE results obtained in 4| . As is well-known, the electronic structure near the 
Fermi energy of the crystalline graphite markedly depends on the weak interlayer interaction 
(see, e.g., [7, ^, 0, Q). Accordingly, the FEC in this case should be sensitive to the specific 
electronic structure. In order to show this, we consider the simplest possible modification 
of graphite (hypothetical AAA stacking) where the three-dimensional energy spectrum was 
calculated analytically in [h]]. Two possible orientations of the applied electric field (along 
and normal to the graphite layers) are of our interest. 



II. FEC OF OPENED CARBON SWNT 



The emitted current density can be written as [3j, [12J 
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Here the field emission is directed along the z-axis, e is the electric charge, h = 2nh the 
Planck constant, e the energy, p momentum, f(e) = [exp(e/A;T) + the Fermi-Dirac 
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FIG. 1: Selected coordinate axes for a rolled (left) and unrolled (right) nanotube. L is the circum- 
ference of the nanotube, l x is a thickness of the graphite layer. 

distribution function, D(e,p x ,p y ) the transmission probability of an electron through a po- 
tential barrier, and v g = de/dp z the group velocity. The integrals are over the first Brillouin 
zone with account of the positivity of v g . Notice that in most cases the using of infinite limits 
in integrals is a good approximation. For parabolic bands v g = p z /m, and this relation is 
widely used in deriving the well-known Fowler-Nordheim equation. In the case of carbon 
nanotubes, two important differences from the generally accepted consideration should be 
taken into account. 

First, an open SWNT has a finite small radius which results in quantization of momentum. 
In this case, the corresponding integrals in ([I]) transform into sums. For example, choosing 
the axes like shown in Fig. [TJone has J f(pi)dpi = ^ f(q)h/k where i — x,y, and l y is the 



circumference of the nanotube, l x the thickness of the graphite layer. For a SWNT there 
exists only a single layer in the x-direction and, accordingly, there is exactly one term in the 
sum for i = x. The number of terms for i — y depends on the tube circumference l y . 

Second, the energy near the Fermi level for a single graphite layer (graphene) is approx- 



move in the yz plane, so that the energy does not depend on p x . Since p y is quantized the 
energy s(p y ,p z ) turns out to be divided into a set of channels with e(hq,p z ) = e q {p z ) where 
q takes integer values. Therefore, the current density ([TJ takes the form 



As is known, the dispersion relation for carbon nanotubes depends on their chirality (see, 



imated by e = ±VF^Jp y + pi, where 




The electrons 




(2) 
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e.g., [III]). For a chiral vector (m,n) it can be written as 



(m — n)/3 + q 

In 



1/2 



(3) 



Generally, there are two symmetric curves with a gap e q g = 2vph((jn — n) /3 + q)/l y . However, 
at certain values of m, n and q the gap turns out to be zero and one gets the linear dispersion 
relation. Therefore, at fixed q there exist one metallic branch and a set of semiconducting 
branches for a SWNT with a given chirality. 

As stated above, the condition v g > imposes restrictions on the limits of integral in (j2J). 
In addition, two approximations will be used. First, we consider the zero-temperature limit 
when the Fermi-Dirac distribution becomes the step function. Second, we suggest that the 



transmission probability is given by the WKB approximation (see, e.g., 
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13] ) in the form 
(4) 



where £ = 87r(2mo) 1//2 /3e/i, y = (eF /Atteq) 1 ' 2 /0, F is the electric field, the work function, 
e the dielectric constant, and we denoted b = exp (— ((p 3 ^ 2 v(y) / F) and d = 3C0 1 ^ 2 t(y)/2F 
for convenience. The functions v(y) and t(y) describe a deviation of the barrier from the 
triangle form due to image effects and can be approximated by [yj] 



v(y) w 1 - y 1 - 



%)« 1 + 0.127?/ 



1.69 



(5) 



Now we are able to calculate the qth term in the sum 
integration in (TSJ) spreads from — oo to 0. One obtains 

4 e F ft 



For the metallic branch the 
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For semiconducting branches the range of integration in (j2J) is (— oo, —e q g /2) and 

3C0 i/2 %^r 



(6) 



J q = Jo exp 



jo exp ( -d-± 



(7) 



2 F 2 

Notice that the dispersion relation enters ($7§ only through the gap. This agrees with the 
well-known fact that the group velocity and the density of states are canceled in the one- 
dimensional case l^ . 

The sum over all branches in ([2]) gives the total FEC. Fig. [2] shows the calculated emis- 
sion current density, which is a current divided by the circumference of a nanotube. For 
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FIG. 2: Current densities vs the diameter of CNTs at the applied field F = 8 x 10 9 V/m. The 
parameter set is 4> = 4.7 eV, v F = 0.83 x 10 6 m/s, a = 2.46 A, ( = 6.83 x 10 9 eV~ 1/2 m 



ly in [4J. Moreover, comparing our results in Fig. [2] with the exact numerical 
4| one can find out a good qualitative agreement. Notice that the quanti- 



-i 



(m, n) SWNT the circumference is defined as l y = a\/m? + mn + n 2 with a being the lattice 
constant. In fact, the main contribution to the sum in ([2]) comes from the first few terms cor- 
responding to branches close to the Fermi level. This is due to the exponential dependence 
of the FEC on the gap. For metallic nanotubes, the leading term is j , so that j™* t ~ l/l y . 
In the case of semiconducting nanotubes, the leading contribution comes from the gth term 
with the smaller gap in ([7]) and, therefore, ~ l/[l y exp(l/l y )]. A similar behavior was 
found numerica 
calculations in 

tative difference is also not great and varies from a few to ten percent depending on the 
diameter of the nanotube. In comparison with [4| our points in Fig. [2] are situated slightly 
lower for metallic nanotubes and slightly higher for semiconducting nanotubes. This differ- 
ence can be explained by at least two reasons. First, we have used the simplified expression 
for the tunneling probability in (jl]) where the image effects were approximated in a standard 
way (see, e.g., Second, as distinct from |4| we consider the zero-temperature limit. 

III. FEC OF CRYSTALLINE GRAPHITE 
A. Noninteracting layers 

In this section, we study a case of noninteracting graphite layers. The layers are oriented 
as shown in Fig. [31 To calculate the FEC we will use the method of independent channels 
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FIG. 3: The location of graphite layers with respect to the electric field. The emission occurs in 
the direction opposite to the electric field. 

described in the previous section. Namely, let us consider the 2D graphite lattice with the 
Born-von Karman boundary conditions applied in the y direction. This gives the natural 
quantization conditions. For layers of a large (infinite) size the sum in ((2]) can be replaced 
by the integral Yliq = (h/h) I J(Py)dPy, so that finally one obtains 



where the relation e g {p y ) = 2vp\p y \ is taken into account. It is interesting to mention that 
this result is very similar to the Fowler- Nor dheim formula 



Indeed, the exponents are exactly the same and the preexponential factors differ only slightly. 
What is important, the independence is equal in both cases. For the interlayer distance 
l x = 3.34 A one can estimate j/j FN = h/ {2vFirmol x ) ~ 0.4. 

B. AAA stacking 

Generally, there are known three possible configurations of crystalline graphite: ABAB... 
stacking sequence of hexagonal layers (Bernal structure), rhombohedral ABC ABC... stack- 
ing, and AAA... stacking when layers of carbon atoms are located directly on top of each 
other [h]] . The AAA stacking is called hypothetical because it has not been observed yet in 
crystalline graphite. However, this configuration is expected in disordered or pregraphitic 
carbon [16]. In this paper, we consider the model of AAA stacking which is the simplest 
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one and allows us to study the effect of interlayer interaction on the emission properties 
analytically. 

Let us consider the interacting graphite layers oriented parallelly to the electric field (see 
Fig. [3]). In the framework of the AAA model the interaction modifies the 3D band structure 
near the Fermi level which can be written in a simple analytical form jlo| 



Here the «o parameter represents the interaction between first-neighboring atoms in a layer, 
a± is related to the interaction between two atoms of the same projection on the yz plane, 
from two neighboring layers, and c is the interlayer spacing. The influence of other parame- 
ters «2 and «3 introduced in [h]] is suggested to be negligible and only linear terms in k • p 
perturbation expansion are taken into account. Actually, the analysis in [lCj shows that the 
maximum effect of the next-to-leading term (breaking the cylindrical symmetry) is of the 
order of five percent. The upper sign in (ITUl) corresponds to the conduction band, and the 
lower sign corresponds to the valence band, aoV3a/2h~ = vp- 

In compliance with ffTUl) Fig. H] represents the Fermi surface of the AAA graphite. The 
Fermi surface is composed of a hole pocket (the valence band, sigh "minus" in (ITU]) ) and two 
half pockets of electrons (sign "plus" in (TTOl ). In our case, the emission occurs along the 
z-axis. Generally, the possible values of the momentum of emitting electrons with respect 
to the Fermi surface can be collected into five different groups. We call them independent 
channels. Solid lines in Fig. H] indicate five possible types of independent channels: (1) for 
an intermediate electron-hole region, (2,3) for holes, and (4,5) for electrons. For finite-size 
layers, quantization of momentum in the xy plane occurs. In this case, the spectrum is 
written as 



where i,j are integer, and p x lies in the region (—7rh/c,7rh/c). The total FEC is a sum of 
all channels. Let us consider these contributions separately. 

1. Hole region 

The hole region is defined as —nh/2c < p l x < irh/2c. As is shown in Fig. H] there are two 
types of channels for the hole region and the channel 1 can be considered as an intermediate 

7 




(10) 




(11) 



electrons 
holes 

electrons" 



FIG. 4: Fermi surface of the simple hexagonal graphite. The central pocket corresponds to the 
hole region, two half pockets correspond to the electron region. The solid lines show five possible 
types of channels for emitting electrons. 

case. Fig. [5] shows all possible one- dimensional dispersion relations for this case. The 
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FIG. 5: One-dimensional dispersion relations for the hole region. 

channel 3 crosses the Fermi surface at two points while the channel 2 does not cross the 
Fermi surface. For each channel = (2e/hl x l y ) J f{e % ^)D{e i ^)de l K The integration here is 
over the occupied states, and is the distance from the Fermi level to the extremum point 
of the branch (see Fig. ED . For the channel 1 A tJ = 0, and the current density is found to be 
equal to jo- Notice that there are only two channels of this type. Analytically, the channel 
2 is denned as \p^\ > rf, where rf = 2aicos(cp l x /fr)/vF- For this channel one obtains 

f=j exp(-dA% (12) 




S 



Here A ?J takes the form A 1 - 7 = vf\p{j\ — 2ol\ cos(cp l x /H), which can be easily found from (TTTT) . 
Replacing the sum ^ j u by the integral one gets 



h/Ac 



A = J d P* J Joexp[-dA(p x ,p y )}dp y = h2 2 d f cvp (13) 

—h/Ac V 

For the channel 3 one has \py\ < rf and j 1 ^ = jo- Like before, one obtains 

h/Ac rj 



2l x l y f i f ■ i 8«ie6 



H = ^ J dp * J jodp » = (14) 

-h/Ac 

2. Electron region 

The electron region is defined as —nh/c < p % x < —nh/2c and Trh/2c < p l x < rrh/c or, 
taken into account a periodicity of the Brillouin zone, 7rh/2c < p l x < 3irh/2c. There are two 
kinds of channels in the electron region with spectra shown in Fig. [6j The channel 4 does 
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FIG. 6: One-dimensional dispersion relations for the electron region 

not cross the Fermi surface and, therefore, the current density turns out to be zero. For the 
channel 5 one has \py\ < \t] l \. This channel contains occupied states below the Fermi level 
and, hence, there is nonzero contribution to the FEC. One gets 



^ = jo[l-exp(-dA^)], 
where A y = 2o>i \ cos(cp x /H)\ — VF\Py\ in accordance with ( fTTI) . Finally, 



(15) 
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dp x J jo{l - exp[-dA(p x ,p y )]}dp y 



h/Ac 



2eb[Adai/it + (l (2tfai) - L (2dai) - l)] 

d 2 cvph 2 



(16) 



where L (x) is the modified Struve function, and Iq(x) is the modified Bessel function. 



3. Resulting FEC 



The total current density is found to be 

l/?T + 

h 2 d 2 cvF 
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E. II 2eb[8dai / it + l (2dai) — ~L (2dai)} . . 



i=2 

Notice that jl ot reduces to (El) for ol\ — 



C. AAA stacking: perpendicular electric field 

Let us consider interacting graphite layers placed normally to the electric field (see Fig. [3]). 
This situation differs markedly from the previous case. Let us denote p = \Jv\ + P'i i n t ne 
dispersion relation in (TlOl) . Quantization of momentum results in the replacement p — > p l K 
There are only two types of channels in this case (see Fig. [7]). As before, let us consider 
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FIG. 7: Fermi surface of the simple hexagonal graphite. The solid lines show two possible types of 
channels for electrons emitting in the rc-direction. 

them separately. 



1. Hole region 

The hole region is defined by —irh/2c < p x < ttH/2c. There are two kinds of channels 
in the hole region with spectra shown in Fig. [81 The channel 1 is defined by p u < 2a\jvp. 
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FIG. 8: One-dimensional dispersion relations for the hole region. 



The current density reads 



fi = j [l-eM-dA ij )\, 



where A 1 - 7 = vpp %3 . One obtains 



2ai/vp 



li=^2~ I io{l - exp[-dA(p)]}2?rpdp 



ineb ( (2a X (i) 



{dhfv\ \ 2 
The channel 2 is defined by p l i > 2a\/vF, and 



1 + exp(— 2aid)(2aid + 1 



(18) 



(19) 



f = j (exp(2aid) - 1) exp(-dA«). 



(20) 



Finally, 



2a\/vF 



I l z f 

h = "fr / jo(exp(2aid) - l)exp[-dA(p)]27rpdp 
o 

4-7re&(2a 1 <i + 1) 
= (dh) 3 v 2 exp(-2aid)(exp(2aid) - 1). 



(21) 



5. Electron region 

In the electron region 7ch/2c < p x < 3irh/2c. There are also two kinds of channels in this 
region with spectra shown in Fig. The channel 3 is equivalent to the channel 5 in the 
previous section. It is defined by p u < 2a\/vp. The current is 



j« = jo[l-exp(-dAtf)]. 



(22) 
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FIG. 9: One-dimensional dispersion relations for the electron region. 



Here A' lj = 2ai — v F p lj (see Fig. [9] and ffTTj) . In this case, 

2a±/vp 

Ji = lj ^ j J {l-exp[-dA(p)]}2npdp 



o 



Aixeb ( (2a 1 d) 



(dh) 3 v 2 F V 2 



— 2a\d + exp(2«i(i) — 1 ) . 



(23) 
(24) 



The channel 4 is defined by p 1 ^ > 2a\jvp. Since it does not cross the Fermi surface the 
current density is equal to zero. 



3. Resulting FEC 



The total current density is the sum of all channels 



IV. DISCUSSION 



Fig. [10] shows the calculated current densities as functions of the applied electric field. 
For comparison, the Fowler-Nordheim curve is drawn. The most important difference comes 
from preexponential factors. As is known, the Fowler-Nordheim theory gives j FN ~ F 2 b at 
all F . In accordance with (|17p the preexponential factor has a different field-dependence. 
At small F one obtains j" ~ Fb. When F increases (which means 2ot\d — > 0) the current 
density comes to j" ~ F 2 b and, finally, j^/j FN 0.4 as was shown in section [III Al Indeed, 



at large x the difference 1q(x) — Lq(x) tends to 2/{irx) while at small x it tends to 1 15]. 



This is clearly seen in Fig. [TT1 where the comparative curves are demonstrated. For j 1 - we 
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FIG. 10: The current densities vs electric field for graphite layers placed parallelly (j") and normally 
(j -1 -) to the electric field with a\ = 0.4 eV, = 5 eV. For comparison the Fawler-Nordheim curve 
(J FN ) is shown. 

have a similar behavior. According to fl25|) . j ± ~ a\Fb at small F, and j 1 - ~ a\F 2 b at large 
F, so that j ± /j FN tends to a constant with increasing F. One can conclude that the bigger 
is the electric field the lesser is the role of the interlayer interaction. The anisotropy of the 
emission from the 3D graphite is also shown in Fig. [TTJ As is seen, j"/i~ L ~ 2.7, that is 
almost a constant in the considered interval of F. Therefore, we obtain three times increase 
in FEC when graphite layers are oriented in parallel with the electric field. 
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FIG. 11: Comparative curves for the current densities vs the electric field characterizing the role 
of the preexponential factor, j" /j 1 - is almost a constant in the considered interval of F. 



It is interesting to discuss the dependence of FEC from the parameter a± which charac 



terizes the interlayer interaction. It was found in 



16] that this parameter is very sensitive 
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FIG. 12: Reduced current density vs the relative interlayer distance x = (c — c*)/c, c = 3.34L4 for 
two field orientations, is markedly more sensitive to x than j'H. 

to the interlayer distance. Based on their results one can approximate 

ai = 18a; 2 - 0.4, (26) 

where x = (c — c*)/c and a\ is measured in eV. As is seen from Fig. [121 there is a strong 
dependence of the FEC from the interlayer distance. The less is this distance the more is 
the emission current. This is valid for both orientations. It would be interesting to check 
this finding in experiments with graphite crystals under pressure. Notice that this result 
follows from the fact that the DOS at the Fermi level (which is of the most importance in 

). As is seen from Fig., [T2], j is more 
sensitive to x than p. Moreover, for oi\ — > one has j 1 — > 0, which follows from the fact 
that the movement of electrons between layers is suppressed in the absence of the interlayer 
interaction. 

Conclusion 

In conclusion, we have found that the band structure of the 3D graphite has a marked im- 
pact on the field emission current. Experimentally, the field emission from carbon materials 
was studied in ^j. Unfortunately, the polycrystalline carbon films used in experiment can 
not be properly described in the framework of our approach because for this purpose we have 
to consider a mixture between different crystalline structures. In fact, the Fermi surface of 
the ABAB structure of graphite is found to be much more complex and, in particular, it does 
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not possess the cylindrical symmetry (see, e.g., [9|, [l6|). In this case, our approach should 
be markedly modified. Besides, many additional factors like the presence of a diamond-like 
phase on the surface of samples and the absence of any information about the local elec- 
tric fields do not allow us to clarify the role of the band structure in this case. Therefore, 
specific emission experiments with graphite single crystals at different orientations of the 
electric field would be of evident interest. 

This work has been supported by the Russian Foundation for Basic Research under grant 
No. 05-02-17721. 
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